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If the dark energy is modelled as a relativistic elastic solid then the standard CDM and ACDM 
models, as well as lattice configurations of cosmic strings or domain walls, are points in the two- 
dimensional parameter space (w,Cg). We present a detailed analysis of the best fitting cosmological 
parameters in this model using data from a range of observations. We find that the y 2 is improved 
by ~ 10 by including the two parameters and that the w = — 1 ACDM model is only the best fit to 
the data when a large number of different datasets are included. Using CMB observations alone we 
find that w = —0.38 ± 0.16 and with the addition of Large-Scale Structure data w = —0.62 ± 0.15 
and log Cs = —0.77±0.28. We conclude that the models based on topological defects provide a good 
fit to the current data, although ACDM cannot be ruled out. 

PACS numbers: 


I. INTRODUCTION 


A Standard Cosmological Model has emerged over the past few years which provides an acceptable fit to all the 
available data 00 It is often known as the ACDM model since the two dominate energy density components of the 
universe at the present epoch are provided by cold dark matter (CDM) and a cosmological constant, A; their densities 
relative to the critical density are denoted fl c and respectively. The model also includes a smaller fraction of 
baryons, given by fib, and density fluctuations created during an epoch of inflation quantified by an amplitude, Ag, 
and spectral index, ns- The expansion rate of the universe at the present epoch is given by Hq = 100ft krnsec^ 1 Mpc -1 
and the optical depth to the epoch of reionization is given by tr. The age of the Universe, to, and the redshift of 
reionization, z Te , can be computed from these parameters. 

Unfortunately, the value of A required to make this model viable is very much in excess of that which is predicted by 
standard arguments from particle theory and much, largely unsuccessful, effort as been expended to try and explain 
the apparent discrepancy between theory and observation. The important feature of A in the observational context 
is that it is constant in space and time and therefore does not cluster into virialized objects, and that it violates the 
weak energy condition, 1 + 3 w > 0 if the pressure and density are related P = wp , which leads accelerated expansion 
in the late universe. The generalized concept of dark energy has been introduced to represent any energy component 
for which 1 + 3w < 0 and whose Jeans length is sufficiently large to prevent clustering. 

A dynamical cosmological constant, often dubbed Quintessence, based on scalar fields has been developed p] as a 
dynamical alternative to A and is the most popular dark energy candidate. If a scalar field, <f >, rolls slowly down a 
sufficiently shallow potential then a fluid with w < —1/3 and a large Jeans length can be arranged by a suitable choice 
of potential. However, in order for the dynamics of the scalar field to be super-damped, the Compton wavelength of 
the field needs to be of order the Hubble radius today which implies that the mass must be ~ H 0 ~ 10~ 33 eV. Such 
small values of m$ are difficult to achieve within particle physics models due to the effects of quantum corrections 0 
and such light particles would naturally couple to the fields of the Standard Model of Particle Physics @. While 
it provides an interesting phenomenological model for dark energy, in many ways Quintessence provides no better 
micro-physical motivation than A. 

The dark energy cannot be a perfect fluid since the sound speed c s is given by 


c 


2 

s 



(1) 


which is < 0 leading to instabilities on the smallest scales. A suggestion :6] has been made that the inclusion of 
rigidity, p, in the fluid, making it an elastic solid, could provide the necessary stability. By definition a solid not only 
supports longitudinal (scalar) perturbations, but also transverse (vector) perturbations. The sound speed of these 
vector modes is given by = p/(p + P) and it can be shown that [j, Q 


dP 4 , 4 9 

= ^ + r*= w+ r- 


( 2 ) 


This equation implies that if p/p is sufficiently large then c 2 > 0 even if w < —1/3. Apart from their role in stabilizing 
the scalar modes, we will not consider the possible effects of these vector modes here, leaving their detailed exploration 
to a future publication :9]. 
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FIG. 1: Allowed parameter space of SDE models. The lines w — 0, w = —1/3, w = —2/3 and w = —1 correspond to the dark 
energy models which scale like CDM, static cosmic strings, static domain walls and A respectively. The point (ui,c 2 ) = (0,0) 
corresponds to CDM, ( — 1/3, 1/5) to cosmic strings with p/p = 4/15 and (—2/3, 2/5) to domain walls with the same value of 
p/p- 


The solid dark matter/energy (SDE) model was first motivated by the desire to model the effects of a lattice 
configuration of non-Abelian cosmic strings CM! , or domain walls [lj, HH which would behave like an elastic solid 
on the scales relevant to cosmological perturbations. On dimensional grounds, the density of static cosmic strings 
scales as p s t T oc a -2 and that for domain walls is given by pp w oc a -1 . Hence, they are natural candidates for dark 
energy with w = —1/3 and w = —2/3 respectively. Recently, it has been shown 0 ] that under some reasonable, 
but not totally general, circumstances p/p = 4/15 for static Dirac-Nambu-Goto strings and domain walls providing 
sufficient rigidity to achieve stability in both cases. Hence, (w, c 2 ) = (—1/3,1/5) represents the cosmic string model 
and (—2/3,2/5) domain walls. What is not totally obvious is that the model presented in ref. [§j, for the which the 
scalar sector of perturbations is completely specified by (w, c 2 ), also naturally models a A fluid with w = — 1 and 
CDM with (0,0). This means that the inclusion of these two parameters provides a parameterization of constant 
dark energy models with four interesting special cases as shown in Fig. ^ The aim of this paper is constrain this 
parameterization using a range of cosmological observations. 


II. METHODOLOGY 

An important feature of dark energy models with w ^ —1 is that, in order to make accurate predictions, one has 
to model the perturbations [l5l | . This has been studied in the context of models based on a scalar field by a number 
of authors flTt 'IT. In the synchronous gauge and for constant w , the perturbations in an elastic solid satisfy SB] 

5s = -(1 + w) (^9 S + ^hj - 3.ff(c 2 - w)5 s , (3) 

9 S = -(1 - 3w)H8 s + —^— k 2 5 s - 2w fc 2 n s , (4) 

1 + w 3(1 +w) 

where H = a/a, h is the trace of the spatial metric perturbation and the anisotropic stress is given by 

n s = ^-(c 2 - w) [-5 S + 3(1 + w){ii - ift)] . (5) 

V — Vi is the change in the anisotropic part of the spatial metric component since the time when the perturbations 
were created, which in the case the adiabatic perturbations under consideration here is conformal time r = 0. It is 
clear that these expression revert to those of a perfect fluid when c 2 = w and, in particular, they model a CDM fluid 
(w, c 2 ) = (0,0). We will denote the energy density in the SDE relative to critical as H s and for w = — 1 this will 
correspond to Da- 
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Parameter 

Prior 


(0.005,0.1) 

a :h 2 

(0.01,0.99) 

9a 

(0.5,10) 

tr 

(0.01,0.9) 

ns 

(0.5,1.5) 

log(10 10 A s ) 

(2.7,5.0) 

w 

(-1.33, 0.33) 

log C s 

(-5,0) 


TABLE I: Table of flat priors. The notation (a, b) for a particular parameter gives the lower and upper bounds allowed in the 
fit. Note that the parameters w and c 2 are also bounded by the constraint that 0 < c 2 < 1. w = —4/3 is the lowest value of w 
allowed by this constraint. 


The equations were incorporated into CAMB 0 , which is a modified version of CMBFAST IQ], and this was extensively 
tested. We note that the viability of topological defect models has been considered before in ref. 01 > where it was 
argued that the perturbations in SDE would behave like those in Quintessence and therefore used the perturbation 
model of ref. 0 - Not surprisingly, this is true near w = — 1, but becomes a poor approximation near w = —1/3. 

This was then incorporated into the October 2004 version of C0SM0MC [5IJ in order to create Markov-Chain-Monte- 
Carlo (MCMC) chains which were used to estimate the confidence limits on the cosmological parameters. Up to 8 
parameters were used in the fit and these were chosen to correspond to some of the known degenerate directions in the 
cosmological parameter space. In particular, we used flbh 2 , fl c h 2 , tr, log(10 10 As), ns, w, logc s and $a the acoustic 
scale defined by the ratio of the sound horizon to the angular diameter distance at the redshift of recombination. The 
intrinsic flat priors, listed in Table QH were chosen to be sufficient broad to incorporate the lines of degeneracy known 
to exist within the space of parameters and so that the maxima were far from the edges. 

We considered six different models, denoted I, II, III, IV, V and VI in the subsequent text. I is the standard 
ACDM model with w = —1; II is the full eight parameter fit to the data allowing both w and logc s to vary; III has 
w = —1/3 (fixed) corresponding to string model with logc s allowed to vary; IV is the same as III but with w = —2/3, 
corresponding to domain walls; the final two models, V and VI, use fixed values of (w, c 2 ), (—1/3,1/5) and (—2/3, 2/5) 
respectively, which are those predicted for static Dirac-Nambu-Goto strings and domain walls with n/p = 4/15 fl4l . 
Models I, V, and VI have six parameters, III and IV have seven, and II has eight. 

The following procedure was used to ensure convergence of the MCMC chains: first, a single chain was run for a 
short period (less than a day) and covariance matrix was computed from it. This was used as input for a further five 
chains which were run until the variance between the chains of the largest normalized eigenvalue of the covariance 
matrix was less than 0.1, a process which took a further one to two days. On a small number of occasions, it was 
necessary to run a further five chains using the covariance matrix of the earlier five chains. We typically obtained 
a variance in the normalized eigenvalue which was between 0.01 and 0.05 for each model considered. Such a low 
variance ensures that each of the five chains agree on the degenerate directions within the parameter space. 

We will use three different kinds of data in our fits. Cosmic Microwave Background (CMB) data from the 1st year 
observations of the Wilkinson Microwave Anisotropy Probe (WMAP) [2j|] (likelihood computed using the routine 
provided by the WMAP team (23j) and the three experiments which observe at substantially higher resolution than 
possible using WMAP. These are the Very Small Array (VSA) j24j, the Cosmic Background Imager (CBI) (2l| 
and the ArCminute Bolometer ARray (ACBAR) J26|. The other two main datasets which were used were Large- 
Scale Structure (LSS) information from the 2DF Galaxy Redshift Survey (2DFGRS), using data in the redshift 
range 0.02 < /c/(hMpc -1 ) < 0.15 where non-linear effects are found to be negligible |23], and measurements of the 
brightness of type la Supernovae (SNe). For the SNe analysis, the set of 42 SNe were used from the Supernovae 
Cosmology Project [28J, at redshifts between 0.18 and 0.23, and 18 from the Calan/Toledo Supernoave Survey (29| . 
at redshifts below 0.1. Using these a sample of 54 were selected which form the primary (fit C) dataset as described 
in lib! . Analytic marginalization was then performed over the absolute magnitude by comparing the effective apparent 
magnitude with luminosity distance evaluated at the relevant survey redshift. In section E] we also investigate the 
implications of the Hubble Space Telescope (HST) key project measurement of h = 0.72 ± 0.08 [3(} and that of 
ttbh 2 = 0.020 ± 0.002 which comes from considerations of Big Bang Nucleosynthesis (BBN) l3l| . 

The interpretation of the LSS data in these SDE models needs to be modified. We have already pointed out that 
one of the defining properties of dark energy is a large Jean’s length, Aj, at the present epoch. In these models, 
Aj ~ c s ro, where tq is the conformal time today and therefore if c s is small (or zero) we can expect the SDE to cluster. 
It is clear that it would be sensible to avoid it being clustered on scales < lOMpc, which implies c s > 0.001, since it 
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Model 

Parameter 

I 

II 

III 

IV 

V 

VI 

n b h 2 

0.023 ±0.002 

0.028 ± 0.003 

0.028 ± 0.003 

0.026 ± 0.003 

0.024 ± 0.002 

0.024 ± 0.002 

n c h 2 

0.103 ±0.014 

0.067 ±0.015 

0.066 ±0.014 

0.079 ±0.016 

0.097 ±0.013 

0.093 ±0.016 

9a 

1.043 ±0.006 

1.050 ± 0.006 

1.050 ± 0.006 

1.048 ± 0.007 

1.044 ± 0.006 

1.044 ± 0.006 

tr 

0.17 ±0.11 

0.46 ±0.11 

0.47 ±0.10 

0.37 ±0.13 

0.16 ±0.08 

0.23 ±0.12 

ns 

0.98 ±0.05 

1.13 ±0.07 

1.13 ±0.06 

1.08 ±0.07 

1.00 ±0.05 

1.02 ±0.07 

log(10 10 A s ) 

3.2 ±0.2 

3.7 ±0.2 

3.7 ±0.2 

3.5 ± 0.2 

3.1 ±0.2 

3.3 ±0.2 

w 

-1 

-0.38 ±0.16 

-1/3 

-2/3 

-1/3 

-2/3 

log c s 

- 

< -1.0 

-1.1 

< -0.7 

-0.35 

-0.2 


0.23 ±0.06 

0.30 ±0.10 

0.31 ±0.06 

0.22 ±0.07 

0.46 ± 0.06 

0.28 ±0.08 


0.77 ±0.06 

0.70 ±0.10 

0.69 ± 0.06 

0.78 ±0.07 

0.54 ± 0.06 

0.72 ± 0.08 

to/Gyr 

13.5 ±0.4 

14.1 ±0.7 

14.2 ±0.3 

13.4 ±0.4 

14.5 ±0.3 

13.7 ±0.4 

%re 

16 ±6 

28 ±3 

28 ±3 

25 ±4 

16 ±5 

19 ±6 

h 

0.76 ±0.08 

0.58 ±0.09 

0.55 ± 0.02 

0.72 ± 0.06 

0.52 ±0.01 

0.66 ± 0.06 

—2 log C 

1453 

1443 

1443 

1446 

1453 

1451 

AIC 

1465 

1459 

1457 

1460 

1465 

1463 

BIC 

1496 

1501 

1494 

1497 

1496 

1494 


TABLE II: Marginalised parameter constraints for the models I-VI using the CMB data only (N = 1384). Note the improvement 
in the value of — 21og£ in going from model I to model II and that the values for V and VI have similar values to I. For sure 
these models cannot be ruled out! The values given for logc s are the upper 95% confidence levels in models II, III and IV. AIC 
and BIC refer to the Akaike and Bayesian Inference Criteria, respectively, which are discussed in section ITYI 


would contribute to dynamical estimates of fl m from rich clusters of galaxies which are sub-critical (see, for example, 
ref. H). 

However, in principle there is nothing to prevent it clustering on larger scales and it is desirable to constrain this 
possibility using the data. The 2DFGRS measures the power spectrum of the galaxies which is usually assumed to be 
related to that of the underlying density field via a linear bias factor, commonly believed to be close to one, via the 
relation P g (k) = b 2 P c (k ) where P{k) is the power spectrum of the relevant component. If there is an additional SDE 
component in the universe which is clustered in some way, then it will also effect the power spectrum of the galaxies. 
In order to attempt to take this into account, we will use 

Pg(k) = b 2 (n s P s {k) + H c P c (fc)^ , (6) 

which makes the, untested, assumption that the biasing of the SDE and CDM is in proportion to their density. This 
will play an important role in our results since CMB observations prefer low values of c s , but these are clearly ruled 
out by 2DFGRS. 

The number of WMAP data points fitted to are 1348 from (899 from temperature and 449 from the temperature 
polarization cross-correlation), 16 from VS A, 13 from CBI, 7 from ACBAR, 32 from 2DFGRS and 54 from SNe. 
Although we will present the number of datapoints, N, involved in our fits in order to attempt to assess the necessity 
or otherwise of the parameters which we have added, some caution should be attached to them since in a number 
of cases the errorbars are correlated and these correlations need to be taken into account. The number of degrees of 
freedom computed in this way is a upper bound. 


III. RESULTS 

We will present constraints on the various parameters in the models I-VI for four different combinations of datasets. 
Firstly, we will just use the CMB dataset (WMAP +VSA+CBI-I-ACBAR) alone, then we will include either 2DFGRS 
or SNe and finally we consider the combination of all three. The detailed results of these different options are presented 
in tables II, III, IV and V. The results of model I can be seen to be compatible with earlier estimates (for example, 
ref. 0) and calibrate the results of model II and models III-VI which are restrictions of it to various special cases. 
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Model 

Parameter 

I 

II 

III 

IV 

V 

VI 

fib ft 2 

0.023 ±0.001 

0.026 ± 0.003 

0.025 ± 0.002 

0.025 ± 0.002 

0.024 ± 0.002 

0.025 ± 0.002 

D c ft 2 

0.115 ±0.007 

0.091 ±0.014 

0.091 ±0.013 

0.095 ± 0.009 

0.097 ±0.011 

0.096 ± 0.009 

9a 

1.042 ±0.005 

1.049 ± 0.007 

1.048 ± 0.007 

1.047 ±0.007 

1.045 ± 0.006 

1.046 ± 0.006 

tr 

0.11 ±0.06 

0.29 ±0.12 

0.27 ±0.11 

0.27 ±0.10 

0.16 ±0.08 

0.23 ±0.10 

ns 

0.96 ±0.02 

1.05 ±0.08 

1.05 ±0.07 

1.04 ± 0.06 

1.00 ±0.04 

1.02 ±0.06 

log(10 10 A s ) 

3.1 ±0.1 

3.4 ±0.2 

3.4 ±0.2 

3.4 ±0.2 

3.1 ±0.2 

3.3 ±0.2 

w 

-1 

-0.62 ±0.15 

-1/3 

-2/3 

-1/3 

-2/3 

log c s 

- 

-0.77 ±0.28 

-0.94 ±0.24 

-0.76 ±0.27 

-0.35 

-0.2 


0.28 ±0.03 

0.28 ±0.06 

0.41 ±0.07 

0.28 ±0.05 

0.46 ± 0.05 

0.28 ±0.05 

D s 

0.72 ±0.03 

0.72 ± 0.06 

0.58 ±0.07 

0.72 ±0.05 

0.54 ± 0.05 

0.72 ±0.05 

to/Gyr 

13.7 ±0.2 

13.6 ±0.4 

14.3 ± 0.3 

13.6 ±0.4 

14.5 ±0.2 

13.7 ±0.4 

■2re 

13 ±4 

22 ±5 

21 ±5 

21 ±5 

16 ±5 

19 ±5 

ft 

0.71 ±0.03 

0.65 ±0.05 

0.53 ± 0.02 

0.66 ± 0.04 

0.52 ±0.01 

0.66 ± 0.04 

—2 log C 

1488 

1477 

1484 

1479 

1490 

1482 

AIC 

1500 

1493 

1495 

1493 

1502 

1494 

BIC 

1532 

1537 

1535 

1530 

1534 

1526 


TABLE III: As for table E] but using the CMB and 2DFGRS datasets (N = 1416). Now the value of logc s is constrained and 
is not just an upper bound. 


ID likelihood curves for w, logc s and ft using model II are presented in Fig. E and these illustrate the main results 
of our analysis. We see that using the CMB data alone the best fitting values of w = —0.38 ±0.16 and that low values 
of c s (< 0.1 at 95% confidence) are preferred. There is an improvement in % 2 = — 21og£ of ~ 10 from model I with 
the inclusion of two extra parameters. This results in a change in the best-fitting value of ft. = 0.58 ± 0.09 which is 
compatible with the HST key project value at the 1 -a level, but is lower than usually assumed. 

There is a reason for these particular values being favoured. As pointed out in ref. 0 the only way the SDE, 
Quintessence and A models can be distinguished using CMB data alone is via the Integrated-Sachs-Wolfe (ISW) 
contribution to the anisotropy, which is created by the decay of gravitational potentials along the line of sight. Since 
the models react to metric perturbations in a different ways, the ISW contribution to the anisotropy is necessarily 
model-dependent. It was shown that SDE models with w = —1/3 and c s = 0 have anisotropy spectra which turn 
down on large angular scales; a known feature of the WMAP (and COBE) datasets. It is for this reason that this 
particular choice of parameters is preferred in this case. This apparent anomaly may be an artifact of, for example, 
imperfect foreground subtraction, but taking the WMAP data at face value one is lead to the conclusion that this 
model is the best-fit. We note that this is somewhat contrary to the standard dogma since the anisotropy spectrum 
in the ACDM model increases on the largest scales. 

The inclusion of the 2DFGRS dataset into the fit shifts the preferred value tow= —0.62 ± 0.15 and dramatically 
constrains logc s = —0.77 ± 0.28. The value of ft = 0.65 ± 0.05 and y 2 is still ~ 10 smaller than in the w = — 1 fixed 
case. The reason for these changes in the preferred values is due to two effects. The first is that the 2DFGRS survey 
measures the cosmological parameter combination f2 m ft ~ 0.2 modifying the preferred values of f l m and ft with the 
consequent effect on w via the relation for 9 a, and also the fact that if c s is low then the SDE contributes to P g (fc) 
in the relevent range of k. A trade-off between this effect, and the desire to keep c s as low as possible in order to 
improve the fit to the large angular scale CMB, results in the preferred value of c s . 

The SNe dataset is known to prefer values of w which are closer to -1 and this is borne out by its inclusion in 
the fit, both in the case when it is only considered in conjunction with the CMB data and also with 2DFGRS. We 
find that w = —0.68 ± 0.14 for CMB ± SNe and w = —0.80 ± 0.15 for CMB ± 2DFGRS + SNe. In these cases the 
improvement in the value of y 2 , when compared to w = —1 fixed, is less, ~ 7, than for CMB data alone and CMB ± 
2DFGRS. 

A selection of 2D likelihoods are presented in Fig. [ 3 ] which illustrate the degeneracies between w and the other 
parameters. These bear out our earlier comments. We see that the preferred value of logc s is independent of w , but 
that the preferred values of D rn and ft are degenerate with w using CMB data alone. This degeneracy is broken by 
the inclusion of data from the 2DFGRS and SNe, so much so that when both are included, the preferred value of w 
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FIG. 2: ID likelihood plots for w, logc a and h. The curves are for the four different combinations of data discussed in the 
text: CMB (dash-dot line), CMB + 2DFGRS (dotted line) and CMB + SNe (dashed line) and CMB + 2dF + SNe (solid line). 
Notice that the preferred value of w moves closer to -1 as more datasets are included. Also that low values of c s are excluded 
by the 2DFGRS data. 



FIG. 3: A selection 2D likelihood contours in model II. The four different columns correspond to the different data combinations 
discussed in the text: CMB (left column), CMB + 2DFGRS (middle-left column) and CMB + SNe (middle-right column) and 
CMB + 2DFGRS + SNe (right column). 
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Model 

Parameter 

I 

II 

III 

IV 

V 

VI 

fib h 2 

0.023 ±0.001 

0.027 ±0.003 

0.029 ± 0.002 

0.027 ±0.002 

0.025 ± 0.002 

0.026 ± 0.003 

n c h 2 

0.108 ±0.011 

0.073 ±0.017 

0.056 ± 0.010 

0.071 ±0.012 

0.086 ±0.014 

0.080 ±0.015 

9a 

1.042 ±0.005 

1.049 ± 0.007 

1.052 ± 0.006 

1.049 ± 0.006 

1.046 ± 0.006 

1.047 ± 0.007 

tr 

0.13 ±0.07 

0.42 ±0.13 

0.54 ±0.07 

0.43 ± 0.09 

0.22 ± 0.09 

0.31 ±0.12 

ns 

0.96 ±0.04 

1.11 ±0.08 

1.17 ±0.05 

1.11 ±0.06 

1.04 ±0.06 

1.06 ±0.07 

log(10 10 A s ) 

3.1 ±0.1 

3.6 ±0.2 

3.8 ±0.1 

3.6 ±0.2 

3.2 ±0.2 

3.4 ±0.2 

w 

-1 

-0.68 ±0.14 

-1/3 

-2/3 

-1/3 

-2/3 

log C S 

- 

< -0.7 

-1.2 

< -0.7 

-0.35 

-0.2 


0.25 ±0.05 

0.19 ±0.05 

0.27 ±0.04 

0.18 ±0.05 

0.40 ± 0.07 

0.22 ±0.06 

fi. 

0.75 ±0.05 

0.81 ±0.05 

0.73 ± 0.04 

0.82 ±0.05 

0.60 ±0.07 

0.78 ± 0.06 

to/Gyr 

13.6 ±0.3 

13.3 ±0.4 

14.1 ±0.3 

13.2 ±0.4 

14.5 ±0.3 

13.5 ±0.4 

■2re 

14 ±5 

26 ±4 

29 ±1 

26 ±3 

18 ±5 

22 ±5 

h 

0.73 ±0.05 

0.74 ± 0.08 

0.57 ±0.02 

0.75 ± 0.05 

0.53 ± 0.02 

0.71 ±0.06 

—2 log C 

1506 

1500 

1509 

1500 

1524 

1508 

AIC 

1518 

1516 

1523 

1514 

1536 

1520 

BIC 

1550 

1558 

1560 

1551 

1568 

1552 


TABLE IV: As for tabled but using CMB and SNe data (N = 1438). The improvement in the fit in going from I to II is less 
marked when compared to tables II and III, but is still reasonably significant. Model V is ruled out in this case, but VI is not. 



Model 

Parameter 

I 

II 

III 

IV 

V 

VI 

Q b h 2 

0.023 ±0.001 

0.024 ± 0.002 

0.027 ± 0.002 

0.025 ± 0.002 

0.024 ± 0.002 

0.025 ± 0.002 

n c h 2 

0.116 ±0.006 

0.102 ±0.014 

0.082 ±0.11 

0.092 ± 0.008 

0.090 ±0.010 

0.093 ± 0.008 

6 

1.043 ±0.005 

1.045 ± 0.006 

1.051 ± 0.007 

1.048 ± 0.006 

1.046 ± 0.006 

1.048 ± 0.006 

tr 

0.11 ±0.05 

0.21 ±0.12 

0.33 ±0.10 

0.29 ±0.10 

0.19 ±0.09 

0.27 ±0.09 

ns 

0.96 ±0.02 

1.01 ±0.06 

1.09 ±0.06 

1.05 ± 0.06 

1.02 ±0.05 

1.04 ± 0.06 

log(10 10 A s ) 

3.1 ± 0.1 

3.3 ±0.2 

3.5 ±0.2 

3.4 ±0.2 

3.2 ±0.2 

3.4 ±0.2 

w 

-1 

-0.80 ±0.15 

-1/3 

-2/3 

-1/3 

-2/3 

logc s 

- 

-0.87 ±0.42 

-0.91 ±0.19 

-0.71 ±0.26 

-0.35 

-0.2 


0.28 ±0.03 

0.27 ±0.04 

0.37 ±0.06 

0.26 ± 0.04 

0.42 ± 0.05 

0.26 ±0.04 

a 

0.72 ±0.03 

0.73 ± 0.04 

0.63 ± 0.06 

0.74 ± 0.04 

0.58 ± 0.05 

0.74 ± 0.04 

to / Gyr 

13.7 ±0.2 

13.6 ±0.3 

14.2 ±0.3 

13.5 ±0.4 

14.5 ±0.3 

13.5 ±0.3 

Zre 

13 ±4 

18 ±6 

24 ±4 

22 ±4 

17 ± 5 

21 ±4 

h 

0.70 ±0.03 

0.68 ± 0.04 

0.54 ±0.02 

0.68 ± 0.04 

0.52 ± 0.02 

0.67 ±0.04 

—2 log C 

1541 

1535 

1552 

1535 

1561 

1538 

AIC 

1553 

1550 

1566 

1548 

1573 

1550 

BIC 

1585 

1593 

1603 

1586 

1605 

1582 


TABLE V: As for table El but using CMB, 2DFGRS and SNe data (N = 1470). 


appears to be independent of them. We also see that w appears not to be degenerate with fib h 2 for the CMB data 
alone, but when all the datasets are used together a degeneracy appears. We will return to this issue in section ITVl 
The models III and IV model dark energy fluids which scale like static cosmic strings (w = —1/3) and domain walls 
(w = —2/3) respectively without restricting to values of c s corresponding to /z/p = 4/15. We see that for CMB data 
alone that the best-fitting parameters for model III are very close to those of model II, and that x 2 ~ 10 better than 
model I for the inclusion of a single parameter. This parameter, logc s , allows a better fit to the data by reducing the 


























































power at large scales as discussed earlier. Model IV is also a better fit to the data than model I. When the 2DFGRS 
data and SNe are included, by themselves and together, model IV becomes a better fit to the data in keeping with 
the results found for model II. 

We have also explored models V and VI which have w fixed and /.i/p = 4/15. These models have the same number 
of parameters as model I and therefore the values of x 2 can be compared directly. We see that in the case of CMB 
data alone that models I and V fit equally as well, whereas model VI provides a marginally better fit. When 2DFGRS 
is also included the fit of model VI is significantly improved over model I; model V is a less good fit, but only by a 
small amount. SNe data rejects w = — 1/3 and this is borne out by our results. Model V provides a significantly 
poorer fit to the data when SNe is included, either by itself or in conjunction with 2DFGRS, whereas model VI is a 
less good fit than model I for CMB + SNe, but is a better fit for CMB + 2DFGRS + SNe. 

IV. DISCUSSION AND CONCLUSIONS 

Our analysis of the SDE parameter space has turned up some interesting features. We have found that the cosmic 
string model with w = —1/3 and p/p = 4/15 is only excluded when SNe data is included. If c s is allowed to vary then 
the poorer fit to the SNe data can be offset by a much better fit to the large angular scale CMB data. It appears that 
the A based model and the domain wall model appear to fit the data equally well at this stage, dependent on ones 
choice of data. Given the dependence of our results on the choice of data and the possible systematic effects which 
may be hidden within the datasets we feel that the main conclusion is that, of the four models discussed in section^] 
only one (the pure CDM model ( w , c 2 ) = (0,0)) can be ruled out with any great certainty, and that the models based 
on topological defects are compatible with the present data. The domain wall model is probably the most favoured, 
but not with considerable significance. 

The current established approaches to the cosmological parameters (CMB, LSS and SNe) are unlikely to separate 
them in the near future. The resolving power of SNe data may be improved by finding more SNe and covering a 
larger redshift range (see, for example, refs. 01 ), but the CMB angular power spectra and linear LSS power 
spectrum at low redshift have limited room for improvement in constraining w 1 1 2l l35l . WMAP polarization data 
may improve the situation by reducing the uncertainty on tr, but the resolving power of total intensity observations 
has largely been exhausted. It is clear that radical alternatives to the dark energy problem are required and many 
new cosmological tests using upcoming data have been proposed. These include: galaxy cluster surveys traced by 
the Sunyaev-Zel’dovich effect and X-ray emission j3CL l3fl : baryonic oscillations in the matter power spectrum |38j : 
cross-correlation of the CMB with measures of LSS 1391: and weak lensing [Ijj. In some cases these are already 
beginning to have an impact (see, for example, refs. I4lll42j~). 

We should comment on the reasons why some earlier analyses might appear to give constraints on w which are 
somewhat closer to -1 than those we present here. Some of this is due to the fact that perturbations in the SDE 
model evolve differently to scalar field motivated models, although this effect is only significant away from w = — 1. 
The main reason is use of priors on h and fl^h 2 from HST and BBN which may contain systematic errors. This is 
apparent from carefully examining the 2D likelihoods presented in Fig. 0 or noting that the best-fitting values of h 
and fib h 2 presented in tables II, III, IV and V clearly depend on the specific model. In order to reinforce this we 
have performed importance sampling on our MCMC chains. This was done by using the covariance matrix from the 
earlier well-converged chains (normalized variance in the eigenvalues typically between 0.01 and 0.05) to run further 
chains generating approximately ten to twenty thousand samples from which to importance sample. 

The best-fitting values of w are presented in table mi for a range of constraints on h to test the sensitivity to 
possible systematic errors in the measurement of the Hubble constant, with and without a prior on fl^h 2 . We see that 
increasing the central value of the external constraint on h reduces the value of w toward -1. This is a strong effect 
when one consider the CMB data alone, but this degenerate direction is already partially constrained when the other 
datasets are included. The addition of the constraint from BBN has an effect in the same direction. If one uses data 
from CMB+2DFGRS+SNe+HST-|-BBN then one comes to the conclusion that w = —0.89 ±0.11, but the possibility 
of their being no systematic effects in these five datasets seems a little unlikely. We believe that constraints from 
two, possible three different datasets, are much more plausible and in just about every possible combination of the 
datasets, models with w significantly different from -1 are favoured. 

We have also assessed whether the inclusion of extra parameters in models II, III and IV are justified by the 
Akaike ^3 and Bayesian Inference Criteria |44), AIC and BIC respectively, as suggested in ref. [45|. The definitions 
are AIC = — 21og£ + 2k and BIC = — 21og£ + fclogTV, where k is the number of parameters in the model. If one 
computes these statistics for models with different numbers of parameters, constrained by the same data, then the one 
with the lowest value of the statistic is preferred. This can be used to justify the necessity for the extra parameters 
included in the fit, although the fact that there are two different statistics prevents one coming to totally unequivocal 
conclusions! It is clear that the improvement in \ 2 for a BIC is much larger since it requires an improvement of ~ 7 
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CMB 

CMB±2DFGRS 

CMB±SNe 

CMB±2DFGRS±SNe 

h 

NO BBN 

BBN 

NO BBN 

BBN 

NO BBN 

BBN 

NO BBN 

BBN 

0.56 ± 0.08 

-0.35 ±0.12 

-0.47 ±0.18 

-0.58 ±0.14 

-0.67 ±0.16 

-0.60 ±0.14 

-0.76 ±0.12 

-0.76 ±0.15 

-0.83 ±0.12 

0.64 ± 0.08 

-0.43 ±0.14 

-0.61 ±0.20 

-0.63 ±0.16 

-0.76 ±0.17 

-0.64 ±0.14 

-0.79 ±0.12 

-0.78 ±0.15 

-0.86 ±0.12 

0.72 ± 0.08 

-0.53 ±0.15 

-0.74 ±0.18 

-0.69 ±0.18 

-0.85 ±0.17 

-0.67 ±0.13 

-0.82 ±0.11 

-0.80 ±0.15 

-0.89 ±0.11 

0.80 ± 0.08 

-0.63 ±0.16 

-0.84 ±0.15 

-0.75 ± 0.20 

-0.93 ±0.15 

-0.71 ±0.12 

-0.85 ±0.10 

-0.82 ±0.15 

-0.92 ±0.11 


TABLE VI: Constraints on w due external constraints on h and PL^h 2 . We have used different constraints on h to represent 
possible systematic effects in its measurement. The HST key project value is h = 0.72 ± 0.08 and BBN implies that fib h 2 = 


0.020 ±0.002. 


for each parameter added, whereas the AIC only requires an improvement of 2. The results presented in tables II, III, 
IV and V are somewhat inconclusive with AIC and BIC leading to somewhat contradictory conclusions. For example, 
using the CMB ± 2DFGRS data model II does considerably better than model I with respect to AIC, but the opposite 
is true for BIC. However, there is some evidence that model IV which has one extra parameter would be selected over 
model I by both AIC and BIC. Clearly more specialized tests would be required to come to an unequivocal conclusion 
on this issue. 

An issue which we have not addressed in this paper is the possibility of intrinsic anisotropy due to the structure of 
SDE. Our equations of motion for the SDE implicitly assume that the distribution of the SDE is initially isotropic, 
and hence the elasticity tensor is specified by two moduli, the bulk modulus (which is related to w) and the rigidity 
modulus fi. It is likely that SDE models based on cosmic strings and domain walls will have point symmetry which 
will require additional moduli. It is possible that these could lead to a reduction in the large angular scale CMB 
anisotropies and the intrinsic alignment believed to exist between different multipoles l46l . Consideration of these 
effects could further improve the fits of these models to the current data and we are currently exploring this. 
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